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Abstract 

The  cold-gas  dynamic-spray  process  is  analyzed  by  numerical  modeling  of  the  impact  between  a  single  spherical  feed-powder 
particle  and  a  semi-infinite  substrate.  The  numerical  modeling  approach  is  applied  to  the  copper-aluminum  system  to  help 
explain  experimentally  observed  higher  deposition  efficiencies  of  the  copper  deposition  on  aluminum  than  the  ones  associated 
with  the  aluminum  deposition  on  copper.  To  properly  account  for  the  high  strain,  high  strain-rate  deformation  behavior  of  the 
two  materials,  the  appropriate  linear-elastic  rate-dependent,  temperature-dependent,  strain-hardening  materials  constitutive 
models  are  used.  The  results  obtained  indicate  that  the  two  main  factors  contributing  to  the  observed  higher  deposition  efficiency 
in  the  case  of  copper  deposition  on  aluminum  are  larger  particle/substrate  interfacial  area  and  higher  contact  pressures.  Both  of 
these  are  the  result  of  a  larger  kinetic  energy  associated  with  a  heavier  copper  feed-powder  particle. 

The  character  of  the  dominant  particle/substrate  bonding  mechanism  is  also  discussed  in  the  present  paper.  It  is  argued  that  an 
interfacial  instability  which  can  lead  to  the  formation  of  interfacial  roll-ups  and  vortices  can  play  a  significant  role  in  attaining 
the  high  strength  of  interfacial  bonding. 

©  2003  Elsevier  Science  B.V.  All  rights  reserved. 
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1.  Introduction 

Coatings  are  being  increasingly  used  to  obtain  the 
required  surface  and  tribological  properties  in  engi¬ 
neering  components.  Such  coatings  are  often  applied 
using  thermal-spray  processes  such  as  high  velocity 
oxy-fuel,  detonation  gun,  and  plasma  spray.  In  these 
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processes,  the  coating  material  is  heated  to  tempera¬ 
tures  high  enough  to  induce  melting.  Consequently, 
the  high  heat  input  to  a  part  being  coated  accompany¬ 
ing  these  processes  can  be  detrimental,  if  the  material 
of  the  part  degrades  when  subjected  to  high  tempera¬ 
tures.  This  problem  is  generally  avoided  in  the  cold- 
gas  dynamic-spray  process.  The  cold-gas  dynamic- 
spray  process,  often  referred  to  as  simply  “cold 
spray”,  is  a  high-rate  material  deposition  process  in 
which  fine,  solid  powder  particles  (generally  1-50  pm 
in  diameter)  are  accelerated  in  a  supersonic  jet  of 
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Nomenclature 

B  body  force  per  unit  mass 

c  interfacial  perturbation  growth  velocity 

c,  imaginary  part  of  the  interfacial  pertur¬ 

bation  growth  velocity 
D  rate  of  deformation 

E  internal  energy 

EY  Young’s  modulus 

J  Yih’s  function  [19] 

L  characteristic  length  scale 

P  pressure 

Q  mass  energy  density  source 

Re  Reynolds  number 

S  deviatoric  stress 

t  time 

T  temperature 

U  characteristic  velocity 

v  velocity 

W  spin  tensor 

Y  yield  strength 

Greek  letters 
a  wave  vector 

£  total  deviatoric  deformation  rate 

£e  elastic  deviatoric  deformation  rate 

£P  plastic  deviatoric  deformation  rate 

At  time  increment 

r  Griineisen  function 

;/  amplitude  of  an  interfacial 

perturbation 

1  instability  wavelength 

p  viscosity 

v  kinematic  viscosity 

p  density 

a  stress 

t  interfacial  perturbation  growth  time 

Subscripts 

R  reference  curve 

0  initial  value 

T  transpose 

Superscripts 
'  trial  quantity 

time  derivative 


compressed  (carrier)  gas  to  velocities  in  a  range 
between  500  and  1000  m/s.  As  the  solid  particles 
impact  the  target  surface,  they  undergo  plastic  defor¬ 
mation  and  bond  to  the  surface,  rapidly  building  up  a 
layer  of  deposited  material.  Cold  spray  as  a  coating 
technology  was  initially  developed  in  the  mid-1980s  at 
the  Institute  for  Theoretical  and  Applied  Mechanics  of 
the  Siberian  Division  of  the  Russian  Academy  of 
Science  in  Novosibirsk  [1,2].  The  Russian  scientists 
successfully  deposited  a  wide  range  of  pure  metals, 
metallic  alloys,  polymers  and  composites  onto  a  vari¬ 
ety  of  substrate  materials.  In  addition,  they  demon¬ 
strated  that  very  high  coating  deposition  rates  on  the 
order  of  5  m2/min  (~300  ft2/min)  are  attainable  using 
the  cold-spray  process. 

A  simple  schematic  of  a  typical  cold-spray  apparatus 
is  shown  in  Fig.  1 .  Compressed  gas  of  an  inlet  pressure 
on  the  order  of  3  MPa  (500  psi)  is  introduced  and  flows 
through  a  converging/diverging  nozzle  to  attain  a  super¬ 
sonic  velocity.  The  solid  powder  particles  are  metered 
into  the  gas  flow  upstream  of  the  converging  section  of 
the  nozzle  and  are  accelerated  by  the  rapidly  expanding 
gas.  To  achieve  higher  gas  flow  velocities  in  the  nozzle, 
the  compressed  gas  is  often  preheated.  However,  while 
preheat  temperatures  as  high  as  900  K  are  sometimes 
used,  due  to  the  fact  that  the  contact  time  of  spray 
particles  with  the  hot  gas  is  quite  short  and  that  the  gas 
rapidly  cools  as  it  expands  in  the  diverging  section  of  the 
nozzle,  the  temperature  of  the  particles  remains  sub¬ 
stantially  below  the  initial  gas  preheat  temperature  and, 
hence,  below  the  melting  temperature  of  the  feed- 
powder  material. 

Because  the  cold-spray  process  does  not  normally 
involve  the  use  of  a  high-temperature  heat  source,  it 
generally  offers  a  number  of  advantages  over  the 
thermal-spray  material  deposition  technologies. 
Among  these  advantages,  the  most  important  appear 
to  be:  (a)  the  amount  of  heat  delivered  to  the  coated 
part  is  relatively  small  so  that  microstructural  changes 
in  the  substrate  material  are  minimal  or  non-existent; 

(b)  due  to  the  absence  of  in-flight  oxidation  and  other 
chemical  reactions,  thermally-  and  oxygen-sensitive 
depositing  materials  (e.g.  copper  or  titanium)  can  be 
cold  sprayed  without  significant  material  degradation; 

(c)  nanophase,  intermetallic  and  amorphous  materials, 
which  are  not  amenable  to  conventional  thermal-spray 
processes  (due  to  major  degradation  of  the  depositing 
material),  can  be  cold  sprayed;  (d)  grain  growth  and 
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Fig.  1 .  Schematic  of  a  typical  cold-spray  system. 


formation  of  the  embrittling  phases  are  generally 
avoided;  (e)  macro-  and  microsegregations  of  the 
alloying  elements  during  solidification  that  accom¬ 
pany  conventional  thermal-spray  techniques,  and  can 
considerably  compromise  materials  properties  do  not 
develop  during  cold  spraying.  Consequently,  attractive 
feed-powder  properties  are  retained  in  cold-sprayed 
bulk  materials;  (f)  “peening”  effect  of  the  impinging 
solid  particles  can  give  rise  to  potentially  beneficial 
compressive  residual  stresses  in  cold-spray  deposited 
materials  [3]  in  contrast  to  the  highly  detrimental 
tensile  residual  stresses  induced  by  solidification 
shrinkage  accompanying  the  conventional  thermal- 
spray  processes;  and  (g)  cold  spray  of  materials  like 
copper,  solder  and  polymeric  coatings  offers  exciting 
new  possibilities  for  cost-effective  and  environmen- 
tally-friendly  alternatives  to  technologies  such  as 
electroplating,  soldering  and  painting  [4], 

The  actual  mechanism  by  which  the  solid  particles 
deform  and  bond  during  the  cold-spray  process  is  still 
not  well  understood.  The  most  prevailing  theory  for 
cold-spray  bonding  postulates  that,  during  impact,  the 
solid  particles  undergo  plastic  deformation,  disrupt 
thin  (oxide)  surface  films  and,  in  turn,  achieve  intimate 
conformal  contact  with  the  target  surface.  The  inti¬ 
mate  conformal  contact  combined  with  high  contact 
pressures  promotes  bonding.  This  theory  is  supported 
by  a  number  of  experimental  findings  such  as:  (a)  a 
wide  range  of  ductile  (metallic  and  polymeric)  mate¬ 
rials  can  be  successfully  cold-sprayed  while  non-duc¬ 
tile  materials  such  as  ceramics  can  be  deposited  only 
if  they  are  co-cold- sprayed  with  a  ductile  (matrix) 
material;  (b)  the  mean  deposition  particle  velocity 
should  exceed  a  minimum  (material-dependent)  critical 


velocity  to  achieve  deposition  which  suggests  that 
sufficient  kinetic  energy  must  be  available  to  plastically 
deform  the  solid  material  and/or  disrupt  the  surface 
film;  and  (c)  the  particle  kinetic  energy  at  impact  is 
typically  significantly  lower  than  the  energy  required 
to  melt  the  particle  suggesting  that  the  deposition 
mechanism  is  primarily,  or  perhaps  entirely,  a  solid- 
state  process.  The  lack  of  melting  is  directly  confirmed 
through  micrographic  examination  of  the  cold-sprayed 
materials  [2], 

Several  research  groups  [5-8],  have  carried  out 
comprehensive  experimental  studies  to  establish  the 
effect  of  various  process  parameters  such  as:  the  type 
of  the  carrier  gas,  the  gas  preheat  temperature,  the 
average  particle  velocity,  the  particle  impact  angle, 
etc.  on  the  deposition  efficiency  (defined  as  the  mass 
percentage  of  the  feed-powder  particles  deposited). 
The  effect  of  these  parameters  can  be  considered  as 
being  well  established  and  rationalized.  However,  the 
effect  of  the  feed-powder  particle  and  the  substrate 
materials  remains  an  outstanding  issue.  For  instance, 
Gilmore  and  co-workers  [8]  (the  results  are  repro¬ 
duced  in  Fig.  2),  found  that  the  deposition  of  copper 
feed-powder  particles  on  an  aluminum  substrate  is 
associated  with  a  significantly  lower  deposition- 
threshold  particle  velocity  and  a  considerably  higher 
value  of  the  deposition  efficiency  in  comparison  to  the 
deposition  of  aluminum-on-copper  at  a  given  average 
particle  velocity.  It  should  be  noted  that  the  data  shown 
in  Fig.  2  pertain  to  the  initial  deposition  efficiencies, 
when  copper  particles  are  impacting  the  aluminum 
substrate  in  the  case  of  copper  deposition  on  alumi¬ 
num  and  when  aluminum  particles  are  impacting  the 
copper  substrate  in  the  case  of  aluminum  deposition 
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Fig.  2.  Deposition  efficiency  as  a  function  of  the  mean  particle  velocity  for  cold  spraying  of  copper-on-aluminum  and  cold  spraying  of 
aluminum-on-copper  [8]. 


on  copper.  Once,  the  substrate  is  covered  with  a  one- 
particle  thick  layer  of  the  coating  material,  further 
deposition  takes  place  by  the  collisions  of  copper 
particles  with  the  copper  substrate  and  by  the  colli¬ 
sions  of  aluminum  particles  with  the  aluminum  sub¬ 
strate.  This  later  stage  of  the  cold-spray  deposition 
process  will  not  be  considered  in  the  present  work. 

The  objective  of  the  present  work  is  to  carry  out 
numerical  simulations  of  the  impact  process  of  a  single 
spherical  feed-powder  particle  with  a  semi-infinite  flat 
substrate  in  order  to  help  elucidate  the  cold-spray 
particle/substrate  bonding  mechanism  and  the  role  of 
the  particle  and  substrate  materials.  Toward  that  end, 
the  computational  method  previously  used  by  Dykhui- 
zen  et  al.  [9]  is  extended  to  study  cold-spray  deposition 
of  copper-on-aluminum  and  aluminum-on-copper. 

The  organization  of  the  paper  is  as  follows:  a  brief 
overview  of  the  equations  governing  the  impact  of  a 
single  feed-powder  particle  with  a  semi-infinite  sub¬ 
strate  is  presented  in  Sections  2. 1-2.3.  Numerical 
solution  of  these  equations  using  a  finite-difference 
based  method  is  discussed  in  Section  2.4.  Details  of  the 
physical  model  used  to  study  the  impact  of  a  particle 


with  the  substrate  are  overviewed  in  Section  3.  The 
main  results  obtained  in  the  present  work  are  presented 
and  discussed  in  Section  4.  The  key  conclusions 
resulted  from  the  present  study  are  summarized  in 
Section  4. 

2.  Computational  procedure 

All  the  simulations  of  the  impact  process  of  a  single 
feed-powder  particle  with  the  substrate  are  carried  out 
using  the  CTH  computer  code  which  was  developed  at 
the  Sandia  National  Laboratory  for  modeling  a  wide 
range  of  solid  dynamics  problems  involving  shock 
wave  propagation  and  material  motion  [10].  In  this 
section,  a  brief  description  is  given  of  the  basic  gov¬ 
erning  equations  implemented  in  the  CTH  code  and  the 
numerical  method  used  for  obtaining  the  solution. 

2.1.  Conservation  equations 

Within  the  CTH  code,  the  total  stress  tensor  is 
decomposed  into  a  hydrostatic  (pressure),  P,  and  a 
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deviatoric  part,  S.  The  evolution  of  the  pressure,  P, 
within  each  material,  during  an  impact  event  is  gover¬ 
ned  by  an  equation  of  state  (e.g.  P  =  f(p,E),  where 
p  is  the  mass  density,  E  the  mass  internal  energy 
density,  and  /  is  used  to  denote  a  function).  The 
deviation  portion  of  the  stress,  S,  within  each  material, 
is  defined  by  a  kinematic  equation  which  combines  a 
yield  criterion,  a  flow  rule  and  a  material  constitutive 
model. 

Within  the  CTH  code,  the  density,  the  velocity  u, 
and  the  internal  energy  density  are  defined  by  the 
following  conservation  equations: 

•  a  mass  conservation  equation 


isotherm  and  Cy  =  ( dE/dT)v  is  the  constant-volume 
specific  heat.  These  equations  are  most  suitable  at 
lower  temperature,  near  the  Hugoniot  and  do  not 
account  for  phase  transitions  such  as  melting,  vapor¬ 
ization,  or  chemical  reactions  and  electronic  excita¬ 
tion.  The  Hugoniot  is  represented  by  a  quadratic 
relationship  between  the  shock  velocity  us  and  the 
particle  velocity,  up  as: 

Up 

us  =  Cs  +  Si  Up  +  S2  —  (6) 

<-s 

where  Cs  is  the  speed  of  sound  while  ,S'i  and  S2  are 
material  constants. 


dp 

1ft 


-pV • b , 


•  three  momentum  conservation  equations 

P^  +  v  •  0  +  IP)  =  - pB , 

•  and  an  energy  conservation  equation 


9(pE) 

dt 


(er  +  IP)  ■  Vv  +  Q 


(1) 

(2) 

(3) 


These  equations  are  solved  in  conjunction  with  an 
appropriate  equation  of  state  and  a  kinematics  equa¬ 
tion.  In  Eqs.  (l)-(3),  t  is  used  to  denote  the  time,  I  the 
second  order  identity  tensor,  B  the  body  force  per  unit 
mass,  and  Q  the  mass  energy  density  source  term, 
while  V  is  used  to  denote  the  gradient  and  V-  the 
divergence  operators  and  a  dot  atop  a  quantity  to 
indicate  a  time  derivative. 


2.2.  Equation  of  state 

The  relationships  between  P,  E,  p,  and  the  tempera¬ 
ture  T  are  defined  by  the  following  Mie-Griineisen 
equations  of  state  [11]: 

P(p,E)=PR(p)  +  pr(E-ER(p)),  (4) 

and 


2.3.  Materials  constitutive  equations 

Two  linear-elastic  rate-dependent  plastic  material 
constitutive  models,  the  Steinberg-Guinan-Lund 
(SGL)  [12,13]  and  the  Zerilli-Armstrong  (ZA)  model 
[14]  both  based  on  a  consideration  of  thermally-acti¬ 
vated  and  drag-controlled  dislocation  dynamics  are 
used  in  the  present  work.  Both  models  take  into  account 
the  effects  of  isotropic  strain-hardening,  thermal  soften¬ 
ing,  strain-rate-dependency  and  pressure-dependency 
of  the  yield  strength,  Y.  In  addition,  the  Steinberg- 
Guinan-Lund  model  accounts  for  the  pressure-  and 
temperature-dependency  of  the  shear  modulus,  G. 

Within  the  Steinberg-Guinan-Lund  model  [12,13], 
the  yield  strength  is  defined  as: 

Y  =  [Tx(ep,  T)  +  T a/(£p)]  (7) 


where  the  thermal  component  of  the  yield  strength,  Yy. 
is  defined  by 


the  function /(ep)  as: 


/(ep)  =  [1  +  /?(£p  +  £/)]"  (9) 

and  the  shear  modulus  as: 


E{PlT)  =  ER{p)  +  Cv{T-TR{p)),  (5) 

where  T  is  the  Griineisen  function  defined  as 
r  =  p~l(dP/dE)p  and  the  subscript  R  is  used  to 
denote  the  quantities  along  a  “reference  curve,” 
usually  the  Hugoniot,  an  isentrope,  or  the  zero-Kelvin 


G(P,  T)  =  Go 


1  +  A—rj-  —  B(T  —  0.02585  eV) 


(10) 


where  YAl  Co,  Ci,  C2,  UK,  YPl  e, , «,  A  and  B  are 
material-dependent  parameters,  ep  the  equivalent 
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plastic  strain,  ep  the  equivalent  plastic  strain  rate  and 
>1  =  p/ p0  the  compression  and  p0  the  initial  density. 
The  temperature  T  is  expressed  in  eV  and  is  defined  as 
kT  ,  where  k  is  the  Boltzmann’s  constant  and  T  is  the 
temperature  in  degrees  of  Kelvin. 

Within  the  Zerilli-Armstrong  model  [14],  the  yield 
strength  is  defined  as: 


tensor  trace.  The  total  deviatoric  strain  rate  e  can  be 
additively  decomposed  into  its  elastic,  £e,  and  plastic, 
£p,  components  as: 

£  =  £e  T  £P  (16) 

The  elastic  deviatoric  strain  rate  is  assumed  to  be 
governed  by  Hooke’s  law  and  can  be  defined  as: 


Y  =  A'  +  (Cl  +  c2v/£p)e'e37’(£p)C47'  +  C5(£pf  (11) 

where  A',  cj,  C2,  C3,  C4,  C5  and  N  are  material-depen- 
dent  parameters. 

When  a  cell  contains  two  or  more  materials,  the 
yield  strength  is  defined  as  the  volume  fraction 
weighted  yield  strength  of  the  individual  materials. 
Consequently,  single  material  cells  containing  voids 
have  a  reduced  yield  strength  since  the  material  volume 
fraction  in  such  cells  is  less  than  unity. 

2.4.  Yield  criterion  and  flow  rule 


Ee  =  ^-(S-WS  +  SW)  (17) 

Z(_7 

where  the  quantity  within  parentheses  is  known  as  the 
Jaumann  co-rotational  stress  rate  and 

W  =  i(Vu-VuT)  (18) 

is  the  spin  tensor. 

The  evolution  of  the  deviatoric  stress,  S,  between 
the  times  t  and  t  +  At,  where  At  is  a  time  increment, 
can  be  generally  expressed  using  a  backward  Euler 
scheme  as: 


The  von  Mises  yield  criterion  is  implemented  into 
the  CTH  code  according  to  which  (plastic)  yielding 
occurs  when  the  magnitude  of  the  deviatoric  stress, 
|S|  =  s/S  ■  S  (a  raised  dot  is  used  to  denote  the  tensor 
inner  product)  reaches  (or  exceeds)  a  material-state, 
temperature  and  loading-rate-dependent  value  of 

\s\  =  yj\y  (12) 

where  the  yield  strength  Tis  given  by  either  Eqs.  (7)  or 
(11).  Upon  yielding  an  associated  flow  rule  is  used  to 
describe  the  subsequent  plastic  deformation.  According 
to  this  flow  rule,  plastic  flow  occurs  in  such  a  way  to 
most  rapidly  accommodate  the  deviatoric  part  of  the 
stress  which  causes  the  plastic  flow  and,  hence,  the 
(deviatoric)  plastic  strain  rate,  ep,  is  defined  as: 


where  I  I  is  used  to  denote  the  magnitude  of  a  tensor. 
The  total  deviatoric  strain  rate,  s,  is  defined  as: 

e  =  £>  -  ±  tr(D)/  (14) 

where  D,  given  as: 

£>  =  1(Vd  +  Vdt)  (15) 


St+At  —  Si  +  St+At  A t  (19) 

where  S  is  defined  by  inverting  Eq.  (17)  as: 

S  =  2 Gse  +  WS  -  SW  (20) 

The  deviatoric  stress  at  the  end  of  a  time  step 
t  =  At ,  St+At,  is  updated  within  the  CTH  code  using 
the  radial  return  procedure  developed  by  Kreig  and 
Kreig  [15],  which  uses  the  total  deviatoric  strain  rate  in 
place  of  the  elastic  deviatoric  strain  rate  to  obtain  a 
trial  deviatoric  stress,  St+ a/,  as: 

S,+At' =  S,' +  2GeAt  (21) 

The  magnitude  of  the  trial  deviatoric  stress  is  next 
relaxed  to  account  for  the  effect  of  plasticity  in 
accordance  with  the  yield  criterion,  Eq.  (12),  as: 

\S,+At\  =  I  Wl  -  2G|sp|  At  =  y[\Yt+ A,  (22) 

Since,  the  shear  modulus  is  a  function  of  T  and  P, 
Eq.  (10),  the  total  strain  rate  a  function  of  v,  Eqs.  (14) 
and  (15),  and  the  yield  strength  is  a  function  of  P,  7] 
|ep|,  and  |ep|,  Eqs.  (7)  and  (11),  Eq.  (22)  can  be  written 
in  a  general  form  as: 

St+At(v,T,P,  |£p|)  =  yf V,+At(P,T ,  |£p|,  |£p|)  (23) 


is  the  rate  of  deformation  tensor,  I  the  identity  tensor,  and  can  be  considered  as  the  governing  equation  for  the 

superscript  T  denotes  the  tensor  transpose,  and  tr  the  magnitude  of  the  (deviatoric)  plastic  strain  rate,  |ep|. 
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Once  |ep  |  is  computed,  the  deviatoric  plastic  strain  rate 
is  determined  using  Eq.  (13),  then  the  deviatoric 
elastic  strain  rate  computed  using  Eq.  (16),  and  finally, 
the  deviatoric  stress  updated  as: 

5,+a,  =  S,'  +  2Gse  At  (24) 

Thus,  Eq.  (24)  in  conjunction  with  Eqs.  (12)— (23)  can 
be  considered  as  the  governing  equations  for  the 
evolution  of  the  deviatoric  stress,  S. 

2.5.  Numerical  solution 

Eqs.  (1)— (5),  (7)  or  (1 1),  (10),  (22)  and  (24),  define  a 
system  of  nine  coupled  partial  differential  and  alge¬ 
braic  equations  with  the  time  t  and  the  spatial  coordi¬ 
nates  x,  y,  and  z  as  independent  variables  and  p,  u,  E,  P, 
T,  Y,  G,  S,  and  £p  as  nine  dependent  variables.  After 
assigning  the  initial  and  boundary  conditions  for  the 
dependent  variables,  the  coupled  system  of  partial  and 
algebraic  equations  is  numerically  solved  within  the 
CTH  code  using  a  finite-difference  based  two-step 
Eulerian  computational  scheme.  In  the  first  step,  the 
Lagrangian  step,  the  computational  mesh  is  attached 
to  the  material  and,  hence,  the  cells  distort  in  order  to 
follow  the  material.  In  the  second  step,  the  Eulerian 
step,  remeshing  is  carried  out  in  order  to  map  the 
distorted  cells  back  to  the  original  Eulerian  mesh. 

3.  Modeling  of  the  particle/substrate 
impact  process 

While  the  angle  at  which  feed-powder  particles 
impact  the  substrate  is  known  to  affect  the  deposition 
efficiency  (e.g.  [8]),  all  the  simulations  carried  out  in 
the  present  work  deal  with  normal-angle  impacts.  This 
was  done  for  the  two  main  reasons:  (a)  normal-angle 
impacts  enable  the  use  of  an  axisymmetric  geometry 
avoiding  the  need  for  computationally  costly  three- 
dimensional  simulations;  and  (b)  the  main  objective  of 
the  present  work  is  to  analyze  the  effects  of  feed- 
powder  particles  and  substrate  materials  and  the  par¬ 
ticle  size,  velocity  and  temperature  on  the  deposition 
process  and  not  the  effect  of  the  angle  of  impact. 

The  substrate  is  modeled  as  a  semi-infinite  200  pm- 
diameter  cylinder,  while  a  20  pm-diameter  sphere  is 
used  to  represent  the  feed-powder  particle.  The  fol¬ 
lowing  two  combinations  of  the  particle/substrate 


materials  are  studied:  (a)  copper/aluminum;  and  (b) 
aluminum/copper.  These  two  combinations  of  the 
particle/substrate  materials  are  considered,  in  order 
to  help  explain  the  experimental  results  displayed  in 
Fig.  2  which  show  that  the  cold-spray  deposition 
efficiency  is  substantially  higher  for  copper  deposited 
on  aluminum  than  the  vice  versa. 

As  discussed  in  Section  2.3,  two  types  of  elastic 
rate-dependent  plastic  materials  constitutive  models 
are  considered  in  the  present  work:  (a)  the  Steinberg- 
Guinan-Lund  [12,13]  and  the  Zerilli-Armstrong  [14] 
models.  The  CTH  materials  database  for  the  Stein- 
berg-Guinan-Lund  model  contains  all  the  parameters, 
except  for  the  strain-rate  parameters,  for  copper  and 
1 100-0-aluminum.  The  strain-rate  dependence  para¬ 
meters  for  the  Zerilli-Armstrong  model  are  available 
within  the  CTH  materials  database.  Unfortunately, 
currently  this  database  contains  the  Zerilli-Armstrong 
model  parameters  only  for  copper.  Hence,  in  order  to 
assess  the  effect  of  strain-rate-dependency  of  yield 
strength  on  the  particle/substrate  impact,  the  following 
particle/substrate  materials  models  have  been  consid¬ 
ered:  (a)  SGL-copper/SGL-aluminum  (no  effect  of 
strain  rate  accounted  for);  (b)  SGL-aluminum/SGL- 
copper  (no  effect  of  strain  rate  account  for);  (c)  ZA- 
copper/SGL-aluminum  (the  strain-rate  effect  account 
for  only  in  copper  particles);  and  (d)  SGL-aluminum/ 
ZA-copper  (the  effect  of  stain  rate  accounted  for  only 
in  the  copper  substrate).  Only  the  results  pertaining  to 
the  cases  (c)  and  (d)  are  presented  in  this  paper  since 
they  are  qualitatively  quite  similar  to  and  quantita¬ 
tively  within  a  few  percent  of  their  counterparts  for  the 
cases  (a)  and  (b).  The  SGL  and  the  ZA  materials 
constitutive  model  parameters  and  the  Mie-Griineisen 
equation  of  state  parameters  for  commercially-pure 
copper  and  1 100-O-aluminum  used  in  the  present 
work  can  not  be  reported  here  since  their  disse¬ 
mination  is  prohibited  by  the  US  Department  of 
State  International  Traffic  and  Arms  Regulations, 
22CFR120-130. 

All  the  calculations  are  carried  out  assuming  that 
the  feed-powder  particle  and  the  substrate  are  initially 
at  room  temperature  while  different  particle  velocities 
within  the  range  400-1000  m/s  are  considered. 

The  computational  domain  containing  a  single 
feed-powder  particle  and  the  substrate  is  discretized 
using  a  mesh  consisting  of  quadrilateral  cells.  To 
ensure  a  fine  cell  size  in  the  particle/substrate  contact 
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Fig.  3.  The  computational  domain  and  the  hierarchical  block 
structure  used  in  the  numerical  simulation  of  the  particle/substrate 
impact. 


region,  an  adaptive  mesh  refinement  strategy  imple¬ 
mented  in  the  CTH  code  is  used.  This  strategy  is 
block-based,  with  the  blocks  (each  containing  the 
same  number  of  cells)  connected  in  a  hierarchal 
manner  with  the  adjacent  ones  and  with  a  maximum 
cell- size  ratio  limited  to  2:1.  An  example  of  the 
hierarchal  block  structure  is  displayed  in  Fig.  3. 
The  interfacial  high-resolution  region  of  the  mesh  is 
created  by  splitting  blocks  midway  along  each  of  the 
coordinate  directions.  This  produces  four  child  blocks 
per  parent  block  for  a  two-dimensional  problem.  Like¬ 
wise,  four  blocks  are  recombined  to  form  a  single 
block  when  such  blocks  are  located  far  away  from  the 
particle/substrate  interface. 

A  symmetry  boundary  condition  is  used  along  the 
vertical  axis  of  the  particle/substrate  assembly,  Fig.  3. 
Sound  speed  based  absorbing  (transmitting)  boundary 
conditions  are  used  along  the  bottom  and  the  side 
surfaces  of  the  substrate.  These  boundary  conditions 
allow  a  mass  flow  into  and  out  of  the  computational 
domain  and  are  generally  used  to  approximate  an 


(a) 


Aluminum 

\  i 


(b) 


jet 

^  Aluminum 


Fig.  4.  The  evolutions  of  shapes  of  a  20  pm-diameter  aluminum  feed-powder  particle  and  a  copper  substrate-crater  for  the  incident  particle 
velocity  of  650  m/s  at  the  times:  (a)  5  ns;  (b)  20  ns;  (c)  35  ns;  and  (d)  50  ns. 
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infinite  or  a  semi-infinite  medium.  Along  the  remain¬ 
ing  portion  of  the  computational  domain,  the  pressure 
is  set  to  zero  within  the  “ghost”  cells  (the  cells  with 
the  void  volume  fraction  equal  to  unity).  The  use  of 
this  boundary  condition  prevents  the  mass  from  enter¬ 
ing  the  computational  domain  while  a  mass  flow  out  of 
the  computational  domain  is  allowed. 

4.  Results  and  discussion 

4.1.  The  evolutions  of  shapes  of  the  particle  and  the 
substrate-crater 

The  simulation  results  of  the  collision  between  a 
single  20  pin-diameter  aluminum  particle  with  the 
incident  velocity  of  650  m/s  and  a  flat  semi-infinite 
copper  substrate  at  four  times  (5,  20,  35  and  50  ns) 
following  the  initial  particle/substrate  contact  are 
shown  in  Fig.  4(a)-(d),  respectively.  According  to 
Fig.  2,  650  m/s  is  slightly  above  the  minimum  critical 
incident  particle  velocity  for  aluminum  cold-spray 


deposition  onto  a  copper  substrate.  The  results  dis¬ 
played  in  Fig.  4(a)-(d)  show  that  as  the  particle/ 
substrate  contact  time  increases,  the  particle  (height- 
to-width)  aspect  ratio  decreases  while  the  substrate- 
crater  depth  and  width  increase.  At  the  same  time,  a  jet 
composed  of  both  the  particle  material  and  the  substrate 
material  is  formed  at  the  particle/substrate  contact 
surface.  When  this  jet  reaches  the  free  surface,  it  forms 
a  lip  which  points  away  from  the  flattened  particle. 
Fracture  of  the  material  at  the  jet  front  as  seen  in 
Fig.  4(b)— (d)  is  controlled  by  a  hydrostatic  tensile  stress 
condition  and,  since  it  is  not  the  focus  of  the  present 
work,  will  not  be  discussed  any  further. 

The  evolutions  of  shapes  of  a  single  20  pm-dia- 
meter  copper  particle  with  an  incident  velocity  of 
650  m/s  and  a  flat  semi-infinite  aluminum  substrate 
at  the  same  four  times  as  in  Fig.  4(a)-(d)  are  shown  in 
Fig.  5(a)-(d),  respectively.  According  to  Fig.  2,  the 
deposition  efficiency  of  copper-on-aluminum  sub¬ 
strate  at  an  incident  particle  velocity  of  650  m/s  is 
~ 5 0-60%\  The  results  displayed  in  Fig.  5(a)-(d) 
show,  like  the  results  displayed  in  Fig.  4(a)-(d),  that 


Fig.  5.  The  evolutions  of  shapes  of  a  20  pm-diameter  copper  feed-powder  particle  and  a  aluminum  substrate-crater  for  the  incident  particle 
velocity  of  650  m/s  at  the  times:  (a)  5  ns;  (b)  20  ns;  (c)  35  ns;  and  (d)  50  ns. 
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as  the  particle/substrate  contact  time  increases,  the 
crater  depth  and  width  increase  while  the  flattened 
particle  aspect  ratio  decreases.  However,  in  this  case, 
the  crater  is  much  deeper  and,  due  to  the  fact  that  the 
particle  is  flattened  to  a  lesser  extent,  the  crater  width 
and  the  final  particle/substrate  contact  surface  area 
(excluding  the  one  associated  with  the  interfacial  jet) 
are  somewhat  smaller.  These  differences  can  be  attrib¬ 
uted  to  the  larger  kinetic  energy  associated  with  the 
heavier  copper  feed-powder  particle  and  to  the  larger 
strength  of  copper  in  comparison  to  that  of  aluminum. 
One  of  the  most  striking  differences  between  the 
corresponding  results  shown  in  Figs.  4(a)-(d)  and 
5(a)-(d)  is  in  the  length  of  the  interfacial  jets.  The 
interfacial  jet  is  much  longer  in  the  case  of  the  copper 
deposition  on  aluminum  than  in  the  case  of  the  alu¬ 
minum  deposition  on  copper. 


Fig.  6.  The  shapes  of  the  feed-powder  particle  and  the  substrate- 
crater  50  ns  following  the  initial  contact  for:  (a)  the  aluminum 
particle  velocity  of  400  m/s;  and  (b)  copper  particle  velocity  of 
500  m/s. 


The  shapes  of  the  feed-powder  particles  and  the 
substrate-craters  50  ns  following  the  initial  contact 
for  the  aluminum  particle  velocity  of  400  m/s  and 
the  copper  particle  velocity  of  500  m/s  are  shown  in 
Fig.  6(a)— (b),  respectively.  These  two  velocities  are 
lower  than  the  corresponding  minimum  critical  par¬ 
ticle  velocity  for  cold-spray  deposition.  By  comparing 
the  results  displayed  in  Fig.  6(a)-(b)  with  the  ones 
displayed  in  Figs.  4(a)-(d)  and  5(a)-(d)  one  can 
observe  that  the  interfacial  jet  is  composed  of  only 
one  material  at  the  sub-critical  particle  velocities.  This 
finding  suggests  that  the  role  of  the  interfacial  jets  in 
removing  the  oxides  and  other  surface  films  from  both 
the  particle  surface  and  the  substrate  surface  is  a 
critical  factor  for  attaining  a  good  particle/substrate 
bonding  (i.e.  a  high  level  of  deposition  efficiency). 

The  time  evolutions  of  the  maximum  temperature, 
the  maximum  plastic  strain  rate,  the  maximum  pres¬ 
sure,  and  the  maximum  von  Mises  equivalent  stress, 
all  at  the  aluminum  particle/copper-substrate  inter¬ 
face  at  the  incident  particle  velocity  of  400,  600,  800 
and  1000  m/s  are  shown  in  Fig.  7(a)-(d),  respectively. 
The  corresponding  results  for  the  copper  cold-spray 
deposition  onto  an  aluminum  substrate  are  shown  in 
Fig.  8(a)-(d). 

The  results  depicted  in  Figs.  7(a)  and  8(a)  show  that 
the  maximum  temperature  at  the  particle/substrate 
interface  briefly  (for  ~  1 1  ns  in  the  copper-on-alumi- 
num  case  and  for  ~20  ns  in  the  copper-on-aluminum 
case)  exceeds  the  melting  point  of  aluminum 
(--940  K)  and  only  at  the  highest  incident  particle 
velocities  (1000  m/s)  used.  The  melting  point  of  cop¬ 
per  (^  1 356  K)  is  never  exceeded.  It  should  be  noted 
that  the  computed  maximum  temperature  values  at  the 
particle/substrate  interface  are  overestimated  due  to 
the  fact  that  the  CTH  code  does  not  enable  modeling 
of  the  transfer  of  heat  from  the  particle/substrate 
interface  region  where  the  heat  is  generated.  A  close 
examination  of  the  temperature  distribution  over  the 
particle/substrate  interface  shows  that  the  conditions 
for  melting  of  the  aluminum  particle  or  substrate  are 
met  over  a  very  small  portion  of  the  interface  (approxi¬ 
mately  14  and  20%  of  the  particle/substrate  interfacial 
area  in  the  aluminum-on-copper  and  the  copper-on- 
aluminum  cases,  respectively).  Also,  the  heat-diffu¬ 
sion  distance  in  the  copper  particle  or  the  substrate 
during  a  typical  collision  time  period  of  50  ns  is  about 
3-4  pm  suggesting  that  heat  conduction  can  play  a 
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Fig.  7.  The  evolutions  of  the  maximum:  (a)  temperature;  (b)  plastic  strain  rate;  (c)  pressure;  and  (d)  von  Mises  equivalent  stress  at  the  particle/ 
substrate  interface  during  the  deposition  of  aluminum-on-copper  at  four  particle  velocities. 


significant  role  in  reducing  the  maximum  interfacial 
temperature  during  the  particle/substrate  collision 
process.  These  findings  suggest  that  melting,  even  if 
it  takes  place,  most  likely  does  not  play  a  key  role  in 
particle/substrate  bonding.  This  is  fully  consistent 
with  numerous  metallographic  examination  results 
which  generally  show  a  cold-worked  grain  structure 
in  the  region  surrounding  the  particle/substrate  inter¬ 
face  in  the  systems  with  good  particle/substrate  bond¬ 
ing  (e.g.  [16]). 


The  results  displayed  in  Figs.  7(b)  and  8(b),  as  well 
as  the  ones  depicted  in  Figs.  4(a)-(d)  and  5(a)-(d), 
show  that,  due  to  a  large  impact  velocity  of  the  feed- 
powder  particle  with  the  substrate,  the  plastic  defor¬ 
mation  is  highly  localized  for  the  region  surrounding 
the  particle/substrate  interface.  Consequently,  a  jet 
is  formed  which  removes  the  oxidized  and/or  other¬ 
wise  contaminated  material  from  the  prior  particle 
and  the  substrate  surfaces.  Consequently,  clean  parti¬ 
cle  and  substrate  materials  are  brought  into  contact 
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Fig.  8.  The  evolutions  of  the  maximum:  (a)  temperature;  (b)  plastic  strain  rate;  (c)  pressure;  and  (d)  von  Mises  equivalent  stress  at  the  particle/ 
substrate  interface  during  the  deposition  of  copper-on-aluminum  at  four  particle  velocities. 


(under  high  pressures)  at  the  particle/substrate  inter¬ 
face  during  the  collision  process.  Formation  of  such 
jets  and  their  role  in  removing  the  contaminated 
material  from  the  particle  and  substrate  surfaces  is 
believed  to  be  one  of  the  major  (necessary)  conditions 
controlling  the  strength  of  particle/substrate  bonding. 

The  results  shown  in  Figs.  7(c)-(d)  and  8(c)-(d), 
show  that,  for  both  aluminum  deposition  on  copper  and 
copper  deposition  on  aluminum,  the  contact  pressure 
and  the  von  Mises  equivalent  stress  values  (the  latter  is 


numerically  equal  to  the  material  yield  strength)  are 
very  high  in  comparison  to  the  yield  strength  in  the 
un-deformed  (kp  =  0)  constituent  materials  at  the 
maximum  interfacial  particle/substrate  interface  tem¬ 
peratures  and  plastic  strain  rates.  For  example,  for  the 
case  of  aluminum  deposition  on  copper  at  the  particle 
velocity  of  800  m/s,  at  the  maximum  interfacial  tem¬ 
perature  (~870  K)  and  at  the  maximum  plastic  strain 
rate  (~660s-1),  zero  plastic  strain  yield  strengths 
of  aluminum  and  copper  are  0.3  x  108  Pa  and 
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0.8  x  108  Pa,  respectively.  The  corresponding  maxi¬ 
mum  interfacial  pressures  and  the  maximum  interfacial 
von  Mises  equivalent  stress  (~4. 1  x  109  Pa  and 
~3.8  x  108  Pa,  respectively)  are,  thus,  higher  than 
the  zero  plastic  strain  yield  strengths  of  the  two  con¬ 
stituent  materials  by  a  factor  between  5  and  130.  This 
finding  suggests  that,  due  to  these  extremely  high 
pressure  and  stress  levels,  it  may  be  perhaps  more 
appropriate  to  treat  the  material  adjacent  to  the  particle/ 
matrix  interface  as  a  viscous  “fluid-like”  material 
rather  than  a  rate-dependent  elastic-plastic  solid  mate¬ 
rial.  Such  a  suggestion  was  previously  made  in  the 
context  of  explosive  welding  (e.g.  [17]),  using  the  same 
argument  as  the  one  presented  above.  By  defining  the 
viscosity,  u,  as  the  proportionality  constant  between  the 
equivalent  shear  stress  and  the  equivalent  shear  rate, 
the  viscosity  of  copper  and  aluminum  under  the  deposi¬ 
tion  conditions  modeled  in  the  present  work  is  esti¬ 
mated,  using  Eqs.  (7)  and  (11),  to  be  between  105  and 
106  P.  Once,  the  material  surrounding  the  particle/ 
substrate  interface  is  assumed  to  be  a  viscous  material, 
one  can  proceed  to  investigate  the  potential  effect  of 
various  fluid-based  phenomena  such  as  interfacial 
instabilities  and  the  resulting  interfacial  waves,  roll¬ 
ups  and  vortices  on  the  strength  of  particle/substrate 
bonding.  The  potential  role  that  such  instabilities 
may  play  in  particle/substrate  bonding  is  discussed 
in  Section  4.3. 

4.2.  Interfacial  bonding  mechanisms 

Many  experimental  and  theoretical/computational 
investigations,  including  the  present  one,  have  docu¬ 
mented  the  formation  of  “surface-scrubbing”  jets, 
high  interfacial  pressures,  large  extents  of  plastic 
deformation  in  the  interfacial  region  and  the  asso¬ 
ciated  increases  in  temperature.  However,  a  generally 
acceptable  explanation  for  the  dominant  mechanism 
for  particle/substrate  bonding  is  still  lacking.  As  dis¬ 
cussed  earlier,  interfacial  melting  may  occur  and 
enhance  the  particle/substrate  interfacial  bonding 
but  it  is  not  believed  to  be  a  major  factor  controlling 
the  bonding  strength.  Formation  of  the  “surface¬ 
scrubbing”  jets  and  high  contact  pressures  are  gen¬ 
erally  considered  as  prerequisites  for  good  particle/ 
substrate  bonding.  However,  these  factors,  while  cri¬ 
tical  for  attaining  clean  particle  and  substrate  surfaces 
and  an  intimate  contact  between  them,  do  not,  per  se, 


offer  any  insight  into  the  character  of  the  dominant 
mechanism  or  mechanisms  of  interfacial  bonding 
under  the  dynamic  cold-spray  conditions. 

The  most  frequently  cited  particle/bonding  mecha¬ 
nisms  under  the  dynamic  cold-spray  conditions  are 
atomic  diffusion,  surface  adhesion  and  plastic  defor¬ 
mation.  In  the  following,  simple  semi-quantitative 
arguments  are  used  to  examine  the  potential  role  of 
these  three  bonding  mechanisms. 

Using  the  Al-Cu  inter-diffusion  coefficient,  Z)CU_A1 
(=10  15  to  10_14m2/s  [18]),  the  atomic  inter-diffu¬ 
sion  distance,  during  the  time  period  of  ~20  ns  over 
which  the  particle/substrate  interface  is  subjected  to 
the  highest  temperatures,  is  found  to  be  only  0.004- 
0.1  nm  (i.e.  a  fraction  of  the  inter-atomic  distance). 
This  finding  suggests  that  the  atomic  diffusion  at  the 
particle/substrate  interface  may  be  excluded  as  a 
dominant  particle/substrate  bonding  mechanisms 
under  the  dynamic  cold-spray  deposition  conditions. 

Adhesion  is  an  atomic  length-scale  phenomenon 
and  its  occurrence  is  controlled  by  the  presence  of 
clean  surfaces  and  high  contact  pressures  to  make  the 
two  surfaces  mutually  conforming.  The  strength  of 
adhesion  then  depends  on  the  (attractive  or  repulsive) 
character  of  atomic  interactions  and  the  crystallo¬ 
graphic  details  of  the  interface.  While  adhesion  can 
be  modeled  using  atomic-scale  simulations,  such  a 
study  is  beyond  the  scope  of  the  present  work.  Never¬ 
theless,  one  can  generally  infer  the  character  of  atomic 
interaction  by  examining  the  corresponding  binary 
phase  diagram.  Specifically,  the  existence  of  a  solu¬ 
bility  limit  is  an  indication  of  repulsive  atomic  inter¬ 
actions,  while  the  formation  of  an  intermediate/ 
intermetallic  compound  is  an  indication  of  attractive 
atomic  interactions.  Complete  mutual  solubility  of 
two  metals  indicates  weak  repulsive  or  attractive 
atomic  interactions.  In  addition  to  the  factors  dis¬ 
cussed  above,  the  overall  strength  of  adhesion-based 
particle/substrate  bonding  depends  on  the  size  of  the 
particle/substrate  interfacial  area.  It  is  generally 
observed  that  stronger  bonding  occurs  when  the 
interfaces  are  wavy  and  when  it  contains  vortex-like 
features  (e.g.  [17]).  This  observation  can  be,  at  least 
partly,  attributed  to  the  associated  increase  in  the 
interfacial  area.  However,  as  discussed  in  the  context 
of  explosive  welding  (e.g.  [17]),  interfacial  per¬ 
turbations  (roll-ups  and  vortices)  can  give  rise  to 
nano/microscale  material  mixing  and  mechanical 
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interlocking,  both  of  which  enhance  the  interfacial 
bonding.  A  more  detailed  discussion  of  these  inter¬ 
facial  instability-based  phenomena  is  presented  in  the 
next  section.  For  now,  the  experimentally  observed 
stronger  bonding  (i.e.  a  higher  deposition  efficiency) 
in  the  case  of  copper  deposition  on  aluminum  than  that 
for  the  aluminum  deposition  on  copper  can  be  tenta¬ 
tively  attributed  to  the  effect  of  higher  contact  pres¬ 
sures  and  a  more  effective  scrubbing  effect  of  the 
interfacial  jet,  both  of  which  increase  the  overall 
strength  of  adhesion-based  interfacial  bonding. 

Plastic  deformation  is  generally  considered  to  be  a 
mechanism  for  energy  dissipation,  i.e.  a  source  of  heat 
in  the  interfacial  region.  However,  as  shown  earlier, 
increased  temperature  in  the  interfacial  region,  is  not 
sufficient  to  enhance  interfacial  melting  or  atomic 
diffusion  to  a  level  at  which  they  act  as  dominant 
interfacial  bonding  mechanisms. 

4.3.  The  role  of  interfacial  instability 

The  discussion  presented  in  the  previous  paragraphs 
suggests  that,  considering  the  very  short  times  during 
which  the  particle/substrate  is  subjected  to  high  tem¬ 
peratures  and  high  pressures,  the  experimentally 
observed  good  interfacial  bonding  at  the  particle 
velocities  exceeding  a  critical  incident  particle  velo¬ 
city,  is  most  likely  the  result  of  some  type  of  nano/ 
microlength-scale  mechanical  material  mixing/inter¬ 
locking  mechanism.  In  the  following,  two  interfacial 
instability-based  mechanisms  are  suggested  by  which 
such  mixing/interlocking  may  occur. 

Interfacial  instability  and  the  resulting  formation  of 
interfacial  roll-ups  and  vortices  can  be,  perhaps,  best 
understood  within  the  context  of  the  Kelvin-Helm- 
holtz  instability  phenomenon.  This  phenomenon  can 
arise  when  two  fluids,  in  contact,  are  moving  at 
different  velocities  in  a  direction  parallel  to  their 
interface.  Instability  can  occur  even  if  the  two  fluids 
have  identical  densities.  When  the  interface  is  sub¬ 
jected  to  a  (non-zero  curvature)  perturbation,  then  as 
one  fluid  flows  around  the  other  a  centrifugal  force  is 
generated.  This,  in  turn,  gives  rise  to  a  change  in 
pressure  which  may  promote  amplification  of  the 
interfacial  perturbation.  As  shown  schematically  in 
Fig.  9,  these  instabilities  may  subsequently  lead  to  the 
formation  of  interfacial  roll-ups  and  vortices.  Such 
interfacial  roll-ups  and  vortices  may  enhance  the 


Material  2 


Fig.  9.  A  schematic  of  the  instability-based  evolution  of  the 
particle/substrate  (Material  1 /Materia  2)  interface  and  the  accom¬ 
panying  formation  of  interfacial  roll-ups  and  vortices. 

overall  strength  of  interfacial  bonding  by  at  least  three 
ways:  (a)  by  significantly  increasing  the  interfacial 
area  available  for  adhesion;  (b)  by  producing  a  fine 
length-scale  mixing  of  the  two  material;  and  (c)  by 
creating  mechanical  interlocking  between  the  two 
materials. 

Interfacial  instabilities  are  typically  understood 
within  the  context  of  inertial  instabilities,  where  flow 
inertia  acts  to  promote  the  instability  and  viscosity  acts 
to  dampen  it.  The  ratio  of  these  two  effects  determines 
if  the  instability  actually  grows  or  decays  and  is 
quantified  by  the  Reynolds  number  as: 


v 
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where  U  is  the  characteristic  velocity,  L  is  the  char¬ 
acteristic  length  scale,  and  v  (=  p/p)  is  the  kinematic 
viscosity.  For  the  situation  considered  here,  a  charac¬ 
teristic  velocity  would  be  the  particle  impact  velocity 
which  is  on  the  order  of  1000  m/s,  and  a  characteristic 
length  scale  would  be  the  particle  diameter  which  is  on 
the  order  of  10  pm.  The  values  of  p  cited  earlier  range 
from  105  to  106  P,  which  translate  to  kinematic  visc¬ 
osities  of  ~100  m2/s  assuming  a  metal  density  on  the 
order  of  1000  kg/m3.  Using  these  numbers,  a  value  of 
Re  ~  1 0  3  is  obtained.  Since  instabilities  arise  when 
the  Reynolds  number  is  substantially  larger  than  one, 
inertial  instability  is  not  expected  to  occur  in  the 
present  situation,  the  rather  rough  estimates  of  length, 
velocity  and  kinematic  viscosity  notwithstanding. 

On  the  other  hand,  the  work  of  Yih  [19]  has 
demonstrated  that  an  interfacial  instability  based  on 
viscosity  differences  between  the  two  flowing  streams 
may  occur,  even  for  vanishing  Reynolds  number.  The 
governing  equations  for  this  viscous  instability  are 
quite  complex  and  their  detailed  analysis  is  beyond  the 
scope  of  the  present  work.  Instead,  a  brief  overview  of 
the  key  equations  is  presented  to  illustrate  the  feasi¬ 
bility  of  this  instability  and  its  potential  role  as  a 
particle/substrate  bonding  mechanism.  Following 
the  analysis  of  Yih  [19],  the  amplitude  of  an  interfacial 
perturbation  )/  is  described  by  the  equation: 

r\  ~  e“(y+ct)  (26) 

where  x!  is  position  along  the  interface,  a  =  2n/X  the 
wave  vector,  and  X  and  c  the  wavelength  and  the 
growth  velocity  of  the  perturbation.  Accordingly, 
the  interface  becomes  unstable  when  the  imaginary 
part  of  the  perturbation  velocity  c,  c\  is  positive  and 
grows  at  a  rate  which  is  proportional  to  the  magnitude 
of  c,.  Yih  [19]  developed  the  following  equation  for  the 
imaginary  part  of  the  perturbation  velocity: 

q  =  ccReJ  (27) 

where  /  is  a  function  of  several  variables  and  is  given 
by  Eq.  (42)  in  reference  [19].  The  magnitude  of  J  is 
generally  of  the  order  of  104  to  105  m2/s.  According  to 
Eq.  (26),  the  amplitude  of  an  unstable  perturbation 
increases  by  a  factor  of  e  over  a  time  period  of 
1/acj.  After  substituting  c;  from  Eq.  (27),  this 
time  period  can  be  expressed  as: 


The  analysis  presented  above  enables  one  to  ascer¬ 
tain  whether  a  perturbation  with  the  wavelength  com¬ 
parable  to  that  of  experimentally  observed  waves 
(ripples)  along  the  particle/substrate  interface  will 
grow  at  a  rate  high  enough  to  form  roll-ups  and 
vortices  during  the  time  scale  of  the  particle/substrate 
collision  event.  For  example,  taking  a  wavelength  of 
1  pm  (smaller  than  the  particle  size  by  more  than  an 
order  of  magnitude),  and  using  the  value  of  the 
Reynolds  number  of  10-4  obtained  earlier  and  a  value 
of  J  =  104  to  105  m2/s  yields  a  time  period  over  which 
an  unstable  perturbation  undergoes  a  significant 
growth  of  t  ~  10-14  to  10-15  s.  Typical  particle/ 
substrate  collision  times  are  on  the  order  of  tens  of 
nanoseconds,  indicating  that  a  1  pm-wavelength  per¬ 
turbation  is  unstable  and  would  grow  during  a  colli¬ 
sion  event.  Furthermore,  Eq.  (28)  reveals  that 
perturbations  with  smaller  wavelengths  (larger  wave 
vectors)  will  grow  even  more  rapidly,  suggesting  that 
this  type  of  instability  can  be  expected  to  result  in  the 
formation  of  interfacial  roll-ups  and  vortices.  Since 
heavier  copper  particles  are  associated  with  larger 
kinetic  energies,  their  collisions  with  the  substrate 
are  expected  to  give  rise  to  formation  of  the  interfacial 
perturbations  with  smaller  wavelengths  in  the  case  of 
copper  deposition  on  aluminum  relative  to  those 
encountered  in  the  aluminum-on-copper  deposition 
case.  This  may  be  one  of  the  reasons  for  the  experi¬ 
mentally  observed  higher  deposition  efficiencies  in  the 
case  of  copper  deposition  on  aluminum. 

The  Yih  interfacial  instability  analysis  [19]  pre¬ 
sented  above  can  further  be  used  to  rationalize  the 
experimentally  observed  stronger  bonding  of  copper- 
on-aluminum  than  that  for  aluminum-on-copper. 
Fig.  2(a)  in  reference  [19]  shows  that  the  magnitude 
of  the  function  J  is  dependent  on  the  viscosity  ratio 
m  =  /q/t2  and  on  the  thickness  ratio  of  the  two  fluid 
streams,  n  =  h\  //i2,  where  subscripts  1  and  2  are  used 
to  denote  the  streams  with  a  lower  and  a  higher 
thickness,  respectively.  For  the  case  of  aluminum 
deposition  on  copper.  Fig.  4(c)  suggests  that  a  thinner 
jet  is  associated  with  the  aluminum  particle,  and  n  has 
been  estimated  as  0. 8-1.0.  Using  the  procedure 
described  in  Section  4.1,  the  aluminum  viscosity  to 
copper  viscosity  ratio  m  has  been  evaluated  as  0.5-1 .0. 
Based  on  Fig.  2(a)  in  reference  [19]  the  corresponding 
magnitude  of  the  J  function  is  assessed  as  (0.5  —  1 .0)  x 
104  m2/s.  For  the  case  of  copper  deposition  on 
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aluminum.  Fig.  5(d)  suggests  that  a  thinner  jet  is 
associated  with  the  copper  particle,  and  n  has  been 
estimated  as  0. 1-0.3.  The  copper  viscosity  to  alumi¬ 
num  viscosity  ratio  m  is  1. 0-2.0.  Based  on  Fig.  2(a)  in 
reference  [19]  the  corresponding  magnitude  of  the  J 
function  is  assessed  as  (1.0— 3.0)  x  105  m2/s.  Accord¬ 
ing  to  Eq.  (28)  and  these  estimates  for  J,  the  growth 
velocity  of  interfacial  perturbations  is  expected  to  be 
larger,  by  at  least  an  order  of  magnitude,  for  the 
copper-on-aluminum  deposition  case  relative  to  that 
for  the  aluminum-on-copper  deposition  case.  Thus, 
the  role  of  interfacial  instability  in  promoting  particle/ 
substrate  bonding  is  expected  to  be  more  pronounced 
in  the  case  of  copper-on-aluminum  deposition  and  is 
another  potential  reason  for  the  observed  higher 
deposition  efficiency.  While  a  detailed  microstructural 
investigation  of  the  differences  in  interfacial  morphol¬ 
ogies  in  the  two  cases  under  consideration  is  currently 
lacking,  a  careful  examination  of  the  scanning  and 
transmission  electron  micrographs  in  references  [3,5], 
seem  to  indicate  that  the  extent  of  interfacial  perturba¬ 
tions  is  somewhat  larger  in  the  case  of  copper  depos¬ 
ited  on  aluminum  than  in  the  case  of  aluminum 
deposition  on  copper.  One  must,  however,  recognize 
that  due  to  high  temperatures  attained  during  the 
particle/substrate  collision  process,  the  microstructure 
in  the  interfacial  region  is  generally  altered  by  recrys¬ 
tallization.  Hence,  experimental  validation  of  the  pro¬ 
posed  interfacial  bonding  mechanism  should  include, 
in  addition  to  a  high-resolution  microstructural  inves¬ 
tigation,  a  chemical  analysis  at  the  same  length  scale. 

Finally,  in  addition  to,  or  instead  of  the  interfacial 
instability  mechanism  described  above,  a  simpler 
topological  bonding  mechanism  is  also  suggested. 
Fig.  5  shows  the  evolution  of  the  particle  and  the 
substrate-crater  shapes  during  an  impact  event  when 
the  particle  is  copper  and  the  substrate  is  composed  of 
aluminum.  A  particularly  important  feature  to  note  in 
this  figure  is  that  the  cavity  formed  by  the  impact  is 
reentrant,  and  that  the  copper  particle  completely 
conforms  to  the  boundaries  of  this  reentrant  cavity. 
This  is  important  because  such  a  cavity,  once  formed, 
prevents  the  copper  material  from  being  separated 
from  the  aluminum  substrate.  For  example,  even  if 
the  two  materials  do  not  adhere  firmly  to  each  other, 
such  a  topological  geometry  would  prevent  separation 
of  the  materials  due  to  the  rivet-like  nature  of  the 
interface.  If  such  a  shape  were  observed  only  when 


deposition  efficiencies  were  greater  than  zero,  and  was 
absent  when  deposition  efficiencies  were  equal  to 
zero,  then  the  presence  of  this  rivet-like  shape  might 
be  the  feature  which  controls  the  bonding  process.  To 
further  explore  this  idea,  the  topologies  of  the  two 
material  systems  under  conditions  where  their  deposi¬ 
tion  efficiencies  are  exactly  zero  are  reviewed.  These 
situations  are  illustrated  in  Fig.  6(b)  for  the  copper-on- 
aluminum  system  and  in  Fig.  4(d)  for  the  aluminum- 
on-copper  system.  In  spite  of  the  different  impact 
velocities,  the  impact  craters  for  these  two  events 
exhibit  remarkably  similar  topologies.  More  to  the 
point,  these  two  craters  appear  to  be  on  the  verge  of 
forming  a  reentrant  shape.  That  is,  the  shapes  are  not 
reentrant,  but  it  seems  clear  that  a  slightly  larger 
impact  energy  would  result  in  a  reentrant  cavity. 
Hence,  it  is  possible  that  the  velocities  at  which  the 
two  systems  explored  in  this  paper  begin  to  show  a 
non-zero  deposition  efficiency  coincide  with  the 
impact  velocities  at  which  reentrant  cavities  begin 
to  form,  creating  a  rivet-like  particle/cavity  assembly. 

5.  Conclusions 

Based  on  the  results  obtained  in  the  present  work, 
the  following  main  conclusions  can  be  drawn: 

1.  The  critical  velocity,  above  which  cold-spray 
deposition  takes  place  is  associated  with  the 
attainment  of  a  condition  for  the  formation  of  a 
particle/substrate  interfacial  jet  composed  of  both 
the  particle  material  and  the  substrate  material. 

2.  Higher  particle/substrate  contact  pressures  and 
better  developed  interfacial  jets,  both  attainable 
through  larger  kinetic  energies  of  the  particle, 
appear  to  be  the  major  factors  controlling  the 
strength  of  interfacial  bonding  and  the  deposition 
efficiency. 

3.  An  interfacial  instability  due  to  differing  viscosities 
and  the  resulting  interfacial  roll-ups  and  vortices 
may  promote  interfacial  bonding  by  increasing  the 
interfacial  area,  giving  rise  to  material  mixing  at  the 
interface  and  by  providing  mechanical  interlocking 
between  the  two  materials. 

4.  A  particle-length  scale  rivet-like  mechanism  may 
also  be  operative  and  its  onset  may  be  linked  with 
the  minimum  critical  particle  velocity. 
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